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We study transport across a line junction lying between two orthogonal topological insulator 
surfaces and a superconductor which can have either s-wave (spin-singlet) or p-wave (spin-triplet) 
pairing symmetry. We present a formalism for studying the effect of a general time-reversal invariant 
barrier at the junction and show that such a barrier can be completely described by three arbitrary 
parameters. We compute the charge and the spin conductance across such a junction and study their 
behaviors as a function of the bias voltage applied across the junction and the three parameters used 
to characterize the barrier. We find that the presence of topological insulators and a superconductor 
leads to both Dirac and Schrodinger-like features in charge and spin conductances. We discuss the 
effect of bound states on the superconducting side of the barrier on the conductance; in particular, 
we show that for triplet p-wave superconductors such a junction may be used to determine the spin 
state of its Cooper pairs. Our study reveals that there is a non-zero spin conductance for some 
particular spin states of the triplet Cooper pairs; this is an effect of the topological insulators which 
break the spin rotation symmetry. Finally, we find an unusual satellite peak (in addition to the 
usual zero bias peak) in the spin conductance for p-wave symmetry of the superconductor order 
parameter. 

PACS numbers: 73.20.-r, 73.40.-c 



I. INTRODUCTION 

Recent theoretical^"^ and experimental^"^ works have 
led to the discovery of a new class of materials called 
topological insulators (TI) . In these materials the surface 
states have a gapless spectrum governed by a massless 
Dirac equation-^°; these states contribute to charge trans- 
port at low temperatures. Such materials could be either 
in two or three spatial dimensions. A two-dimensional 
TI will host a one-dimcnsional gapless edge states while 
three-dimensional (3D) TI will host gapless states pro- 
tected by time-reversal symmetry on its two-dimensional 
surface. The 3D TIs can be classified as strong or weak 
depending on whether the number of Dirac cones is odd 
or even, and this number is determined by a topological 
invariant. The odd number of Dirac cones on the surface 
of a strong TI is protected against time-reversal invariant 
(for instance, non-magnetic) perturbations for topologi- 
cal reasons^'^'^. For materials such as Bi2Te3 and Bi2Se3, 
surfaces have been found which have a single Dirac cone 
near the F point of the 2D surface Brillouin zone^"^. 
Many interesting features of the surface Dirac electrons 
have been studied^'^"'^""'^®. Some of these studies involve 
interfaces created in a TI using proximate magnetic or su- 
perconducting materials or gate voltages^^'-'^'^'^'^'^^"^®'^^. 
Junctions of different surfaces of a TI (in some cases 
separated by a geometrical step or a magnetic domain 
wall)^^"^® or of surfaces of a TI with normal metals or 
magnetic materials^^ have also been studied. However, 
junctions of multiple TI surfaces with a superconductor 
have not been studied in detail so far; this is what we 
aim to do in this paper. 

The problem of electron transmission across a junc- 
tion of a normal metal (NM) and a superconductor (SC) 
has been extensively studied for many years^°"^^. It is 



well-known that the sub-gap transport in such junctions, 
for small barrier strengths, is governed by Andreev re- 
ficction while for large barrier strengths they reflect the 
sub-gap quasiparticle density of states (DOS)^^''^". Since 
the presence of localized quasiparticle edge states below 
the gap depends on the order parameter symmetry, such 
transport measurements, in the strong barrier limit, pro- 
vide us with a tool for determining the pairing symme- 
try of the superconductor. In contrast, in the weak bar- 
rier limit, it was demonstrated that the conductance of a 
NM-SC junction can be more than the NM-NM conduc- 
tance owing to Andreev refiection^^. More recently, hy- 
brid superconductor-semiconductor (Bi2Te3) devices has 
been successfully fabricated and conductance measure- 
ments in such systems has been carried out^^. Wc note 
that these systems are somewhat akin to the junctions 
that we study in this work. However, the junction of a 
TI and a SC is expected to be more complex than its NM- 
SC counterpart. The reason for this complexity arises out 
of the fact that two-dimensional (2D) surface states of a 
TI display spin-momentum locking; hence scattering at 
a junction that changes the electron momentum couples 
different components of the electron spin. In contrast, 
the particle and hole are coupled to each other on the SC 
side. Hence a treatment of the transport between a TI 
and a SC necessarily requires us to use a four-component 
spinor formalism describing both spin and particle-hole 
degrees of freedom^''"^^ . 

In this work, we will analyze the transport properties 
for the system shown in Fig. 1. The system consists of 
two orthogonal TI surfaces, called TI-1 and TI-2. There 
is a two-dimensional SC surface which can be formed ex- 
perimentally by depositing a 2D superconducting film on 
the surface of a 3D insulator. In what follows, we shall 
consider both singlet s-wavc and triplet p-wave pairing 



2 



symmetries for the superconducting film. The two TIs 
and the SC are separated by a line junction. We will 
consider a Dirac electron incident on the line junction 
from TI-1 with arbitrary angle of incidence and study 
its reflection (normal and Andreev) back to TI-1 and its 
transmission (normal and Andreev) into the TI-2 and the 
superconducting film. The main results that we obtain 
from such an analysis are the following. First, we develop 
appropriate boundary conditions for studying transmis- 
sion across such a junction involving Dirac electrons in 
the TI and Schrodinger electrons in the superconductor. 
The general time-reversal invariant boundary condition 
is found to involve three parameters which can be in- 
terpreted as the strengths of three barriers close to the 
junction on the TI-1, TI-2 and SC sides. Second, using 
these boundary conditions, we compute the charge and 
spin conductances of such a junction as functions of the 
barrier strengths and the bias voltage and thus compare 
and contrast the properties of sub-gap transport in these 
junctions with their conventional counterparts. In par- 
ticular, we find that for Xi = X2 = Oi these junctions 
never reach the maximum charge conductance value 2Go 
(where Go is an unit of conductance defined in Eq. (35)) 
found in conventional NM-SC junctions. Third, we find 
that, in contrast to all conventional NM-SC junctions, 
TI-SC junctions can be used to distinguish between the 
different spin states of the Cooper pair of a triplet p-wave 
superconductor. We show that such a property stems 
from the spin-momentum locking of the electrons on the 
TI surfaces. Fourth, we find that for ap-wave SC, the dif- 
ferential conductance has a zero-bias peak similar to the 
conventional NM-SC junctions; however for finite biases, 
the conductance both oscillates and decreases as the bar- 
rier strength increases which is to be contrasted with the 
monotonically decreasing nature of sub-gap conductance 
in conventional NM-SC junctions. Thus these junctions 
display both Dirac-like and Schrodingcr-like characters. 
Finally, we find that there is a non-zero spin conductance 
for one particular component of the spin and two possible 
spin states of the triplet Cooper pairs. Further, we find 
that the differential spin conductance displays an unusual 
satellite peak away from zero bias and study the behavior 
of this peak with chemical potentials of the TI and SC 
surfaces and the barrier strengths of the junction. 

The detailed plan of this paper is as follows. In Sec. II 
we write down the Hamiltonian for electrons on the 
topological insulator surfaces and the superconductors 
for both s— and p— wave pairings and discuss their ba- 
sic properties which will be useful in subsequent analy- 
sis. This is followed by Sec. Ill where we chart out the 
boundary conditions appropriate for transport through 
the junction. We utilize current conservation to dis- 
cuss how the current gets converted from single quasi- 
particles (electrons or holes) near the junction to Cooper 
pairs deep inside the SC. Next, in Sec. IV, we use the 
boundary conditions to compute the transmission ampli- 
tudes of electrons from the TI-1 into the TI-2 and the 
SC and hence the charge conductance of the junction for 



arbitrary bias voltage, chemical potential difference, and 
barrier potential parameters. We discuss the obtained 
results in details in Sec. V. This is followed by Sec. VI 
where we study the spin transport through a p-wave SC 
with different pairing symmetries. Finally, in Sec. VII, 
we present a discussion of our results and some possible 
experiments to test our theory. 
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FIG. 1: Schematic diagram of the proposed system with a line 
junction. Topological insulator surfaces 1 (lying in the x — y 
plane with y < 0) and 2 (lying in the x — z plane with 2 < 0) 
and a superconductor (lying in the x — y plane with y > 0) 
meet at a junction given by the line y — z — 0. We indicate 
several processes which can be triggered by an electron inci- 
dent from the TI-1 side: normal or Andreev reflection back to 
TI-1, normal of Andreev transmission to the TI-2 sides, and 
four possible BdG quasiparticle transmissions on the SC side 
(which eventually decay to zero and get converted to Cooper 
pairs) . 



II. HAMILTONIANS 

The Hamiltonian for the electrons on the surface of a 
TI is given by^^ 

Hn = J -^^'^l^^'^pf^ -a xk- (1) 

where vp is the Fermi velocity, the unit vector fi points 
normal to the 2D surface, fc,;, kj are momenta in the 
2D plane, and ipj: = (V'jr^i ^fcj.)"^ ^ two-component 
spinor. For TI-1, n — z and (i,j) ~ (2;, y), while for 
TI-2, n = y and = {x,z). The energy-momentum 

dispersion and the eigcnstates on the TI-1 side are given 
by 

E = -/ii ± hvi^kl + fc2, 
- (I/V2) [1, T^e'''f, 
e'^^ = ik,+tky)/^kl + kl, (2) 

where vi is the Fermi velocity, and Oi lies in the range 
[0,7r]. In our calculations, we will consider the band 
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which has the + sign in the expression for the energy 
by choosing > and working at energies \E\ < fii. 
On the TI-2 side, the dispersion and the eigenstatcs arc 
given by 



E 



-/i2 + hv2\Jkl + fc2, 

(I/V2) [V1+COS02, Vl-COS02]^, 

(fc^ - ik,)/^kl + kl , 



(3) 



where V2 is the Fermi velocity. To keep the discussion 
simple we choose /ii = = fJ-Ti- Also, our numerical 
results will be obtained for the case vi = V2 although we 
have retained the general form {v2 < wi ) in the important 
analytical expressions. 

The wave functions on the SC side are described by 
a four-component Bogoliubov-de Gennes (BdG) spinor 

whose upper two components correspond to particles 
and lower two components correspond to holes. Namely, 
^' = [ip^, ipi, if}*, '0* ]"^- (Note that many papers use 
the Nambu convention, with = [i/j^, ipi, — i'l, i^* ■ 
In that convention, for example, factors of ia^ will not 
appear in Eq. (5)). In terms of the Pauli matrices a^'V'^ 



and 



which act on the spin and particle-hole com- 



ponents respectively, the Hamiltonian on the SC side can 
be written as 



H. = 



dx I dy ^'^(x, y) 
'0 



^J'Sc T 



2m 

A{x,y) *(.T,y), 



(4) 

where A{x, y) is called the pair potential. For our anal- 
ysis of an s-wave SC (in which Cooper pairs form a spin 
singlet), we will take the pair potential to be of the form 



A{x,y) 







(5) 



A{x,y) 







-Aor{k){d-s*)iay 



\of{k){d ■ a)ia 




(6) 



where d is a unit vector with real components, and f{k) 
is defined below. Physically, d governs the spin pairing 
of the Cooper pairs in the p-wave SC. For d = x, y and 
z, a consideration of the matrix {d ■ a)ia''^ shows that the 
Cooper pairs are in the spin states ft — iij tt + iii and 
+ It, respectively. (In principle, d could depend on 
the momentum k. However, in this work we are going 
to study systems with constant d). In Eq. (6), /(fc) is a 
dimensionless linear function of k defined as 



f{k) = -i{a.j,dx + aydy), 
f*{k) = i{a*^d^ +aldy), 



(7) 



where constants which may be complex. [In 

this paper we will consider wave functions in the SC 
which are of the form e'^i'^^-^+'^yv) where k^ is real but 
ky may be complex. Hence there will not be any simple 
relation between f{k) = a - k and f*{k) = —a* ■ k in gen- 
eral]. In both Eqs. (5) and (6), Aq will be assumed to be 
a real parameter with dimensions of energy. 

In a SC, the nature of the energy-momentum disper- 
sion depends crucially on the symmetry of the supercon- 
ducting pair potential. The dispersion is given by 



E{k) = ± 



h^kl+k^) n2 
-/Use 



2m 



(8) 



for an s-wave SC. For a p-wave SC, we can use the fact 
that (d • a)iay[d ■ (j*)i(jy = 1 to show that 



± 



h^{kl 



2m 



Ag (d- fc)(d* • k). (9) 



The important difference to note between the dispersions 
for these two kinds of pair potentials is that the super- 
conducting gap is always isotropic for s-wave symmetry 
while in the case of p-wave, it is isotropic only when 
(d • fc)(d* • k) is invariant under rotation in x — y plane. 



III. PROBABILITY AND CHARGE CURRENTS 
AND BOUNDARY CONDITIONS 



In this section, we will first discuss the equations of 
motion. Then we will define a probability density pp and 

current Jp, and a charge density pc and current Jc- It is 
well known that the corresponding continuity condition 
for superconductors necessitates conversion between par- 
ticle and condensate currents as discussed for s-wave in 
Rcf. 30. Here we shall discuss the continuity condition 
for both s- and p-wave superconductors using the four 
component formulation. 

Given the two-component spinor ip, let us define = 
-0*. For a TI with a unit normal n, the equations of 
motion for '0 and 4> are given by 

ihdtip = [hvpfi ■ S X (— iV) — P-Ti\ Ip, 

ihdt(t) = [fi-vpn ■ <?* X (— «V) + ^ti] 0- (10) 

For an s-wave SC, the equations of motion are 




Psc j 'P + Aoicr^0, 

- psc\^ - AoiaV-(ll) 
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For a p-wavc SC, the equations of motion are 



ihdtip 



ihdf 



2m 



Aof{k){d-a)i<jycj), 

( fi'v^ \ 

Aor(fc)(^r•a*)^a>. 



(12) 



In deriving the equations of motion for (jj in Eqs. (10- 
12), we have used the fact that 4> is the complex conjugate 
of ip. We will now begin to treat ip and 4> as independent 
variables which are not related by complex conjugation; 
for example, when solving the equations of motion, we 
will assume that ip and (j) depend on the time t in exactly 
the same way, i.e., as e"*^*/''. The upper two and lower 
two components of the BdG spinor 5' will be given by tp 
and (j) respectively. 

The probability and charge densities are given by 

= e(ip^ - 4,^) ^ e*V^*. (13) 

We will now look for currents Jp and Jc which can satisfy 

the continuity equations of continuity, namely, dtPp + V ■ 

Jp = and dtPc + V • Jc = 0. In the TIs, we find that 
the currents 

Jp ~ vpli^'^n X ail; + (jy^h x a*(j)], (14) 
Jc ~ evpi'ip'^n X (j-i/) — (jy'n x (?*</)], (15) 

satisfy the respective continuity equations. 

s-wave SC: In an s-wave SC governed by Eq. (5), we 

find that 



Jp = — /m(V'^VV') - —Im{(t)^V(j}) 
m m 

= — /m(^'tr^V*) 



(16) 



(where Im denotes the imaginary part) satisfies the con- 
tinuity equation with pp. However, it is well known that 



X = —Im{tp''yi)) + —Im{(j)^V(t)) 
m m 

m 



(17) 



does not satisfy the continuity equation with pc- Instead, 
if we define the Cooper pair currcnt'^'^ 



J 



pair 



4eAn /"^ 

y dy' Re[r{x,y')ayct>{x,y% 



(18) 

then the total charge current in the superconductor, de- 
fined as Jtot = Jc + Jpair, satisfics dtPc + ^ ■ Jtot = 0. For 



a state which has definite values of the energy E and the 
momentum k^, we have dtPc = and dxJx.tot = 0; hence 
Jy.tot = Jy.c + Jy,pair wiU also be independent of the y 
coordinate. However Jy^c and Jy^pair will separately vary 
with y; as y goes from to oo, Jy c will go from a finite 
value to while Jy^pair will go from to a finite value. 
In other words, the single electron or hole current Jy ^ 
will gradually get converted to the Cooper pair current 
Jy^pair as y increases 

p-wave SC: In a p-wave SC governed by Eq. (6), we find 
that if f{k) = —i{axdx + o-ydy), then 



J rj 



-/m(^'V^V*) 



2Ao + - 

-I- X — ^ Re[axil}'{d ■ (7)10^(1)] 

2A 

+ y ^ Reiayip^d- a)iay<P] (19) 

satisfies the continuity equation with pp. (Note that the 
expression for Jp in the p-wave case contains a term pro- 
portional to Aq, unlike the expression in the s-wave case 
which docs not have such a term). The charge current 
Jc defined in Eq. (17) again fails to satisfy the continu- 
ity equation with pc- But if we define the Cooper pair 
current 



Jnfj 



2eAo 



x / dy' Re[^Hx,y')id-<j)ayfik)4{x,y') 

^{f{k)i,\x,y')){d- a)oy<l,{x,y')] 



(20) 



the total current Jtot 
0. 



Jc + Jpatr satisfies dtPc+V-Jt, 



Boundary Conditions 

We now turn to the boundary conditions which need to 
be imposed at the junction to ensure that the component 
of the probability current normal to the junction is con- 
served. We will denote the currents on the three sides of 
the junction by Jpi, Jp2, and Jp^. On the TI-1 and TI-2 
sides, the incoming currents are given by (y • Jpi)y_>o- 
and (z • Jp2)z~^o- respectively, while on the SC side, the 
outgoing current is (y • Jp3)j,->o+- The current conserva- 
tion condition is therefore 

{y ■ Jpl)y^O- — {z ■ Jp2)z-^0- = {y ■ Jp3)y^0+- (21) 

Using the expressions for Jp given in Eqs. (14), (16) and 
(19), we find that Eq. (21) implies 



Imi^lT'dy^s) (22) 



TO 
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if the SC is s-wave, and 



IV. CONDUCTANCE CALCULATIONS 



(23) 



if the SC is p-wave. Now onwards, we shall use the 
four-component spinor ^ instead of the two-component 
spinors ip and (j). 

To find the general boundary condition at the junc- 
tion which satisfies Eqs. (22) and (23), let us assume that 
there are three barriers which lie on the TI-1, TI-2 and 
SC sides of the junction and are located very close to the 
junction; we will model them as (5-function barriers with 
dimensionlcss strengths xij X2 and X3 respectively. We 
will assume that these barriers are invariant under time 
reversal, for instance, that they do not involve any mag- 
netic fields. Then an analysis similar to that in Ref. 29 
will give the following boundary conditions at the junc- 
tion for the case of an s-wave SC, 



mvi 



c [Af(xi)*i+/?A/t(x2)«'2], 



2X3*3 



[M(xi)*i-/3Mt(x2)*2 



where M(x) — cosx — i sinxc^ ® 



(24) 



where /3 is related to th e ratio of the Fermi velocities in 
TI-1 and TI-2 as /? = yf^hl. 

For the case of a p-wave SC, a simple generalization of 
the boundary condition in Eqs. (24) which conserves the 
probability current at the junction is given by 



*3 = c[M(xi)*i+/3Mt(x2)*2], 
dy^i - 2x3*3 

Ao_ / -ay{d-G)ay 
hvi I al{d-a*)ay 



mvi 



*3 



= ^a"(»T^ [A/(xi)*i-/3Aft(;^2)*2].(25) 

It can be shown by a simple calculation that the bound- 
ary conditions in Eqs. (24-25) also satisfy charge current 
conservation at the junction. 

Note that Eqs. (24-25) contain a real dimensionlcss pa- 
rameter c. A precise determination of the value of c re- 
quires a microscopic knowledge of the junction and is be- 
yond the scope of the present work. A similar parameter 
appears in the study of junctions in other systems^^'^^. 
In the limits c — > or oo, the SC gets decoupled from the 
TIs and the system reduces to one which only involves 
two TI surfaces. (The problem of two TIs with a junction 
has been studied in Ref. 22). In this paper we will set 
c = 1 in all our numerical calculations. 



We are interested in the charge transport at a sub-gap 
applied voltage between the TIs and the SC. Our calcula- 
tion will proceed as follows. Given an electron incident on 
the junction from the TI-1 side with unit amplitude at an 
angle of incidence 9i and energy E and using the bound- 
ary conditions (Eqs. 24 and 25), we shall compute the 
amplitudes of eight other wave functions, namely, nor- 
mally refiected electrons and Andreev refiected holes on 
the TI-1 side (with amplitudes and va respectively), 
normally transmitted electrons and Andreev transmitted 
holes on the TI-2 side (with amplitudes and i^i). and 
four electron and hole-like BdG quasiparticle wave func- 
tions on the SC side with amplitudes ii, t2, ts and t4. 
Note that in what follows we shall set the zero energy at 
the middle of the superconducting energy gap. 



A. Wave Functions in the TIs 

To account for a possible Andreev reflection and An- 
dreev transmission in TI-1 and TI-2, we have to write 
down the wave functions in the TIs as four-component 
BdG spinors. Due to translational invariance along i, all 
the wave functions will have a factor of 6**^="^. Also, all 
the excitations will be taken to have an energy E which 
means that there will be a common factor of e"*^*/''. 
With r^v and rA being the amplitudes for normal and 
Andreev reflections, we can write the wave function in 
TI-1 as 



1 



*i = ^ 



1 



V2 1 I' + V2 




/ 1 



' 



V 



rA 
'V2 



(26) 



Here, {kx,ky) is the momentum of the incident electron 
on TI-1 which is related to {£,61) as in Eq. (2). The 
normally refiected electron has a momentum {kx,—ky) 
and the Andreev reflected hole has momentum {kx,kyfi) 
where kyh = ^(/^t/ - EY/{hviY - kl and e''^!" = 
{kx + ikyh)/ ^Jk^ + kyi^. For the case E = 0, the hole 
will have a momentum {kx,ky) just like the incident 
electron, but its group velocity Vg = will be op- 

posite to the incident electron's group velocity. This is 
called retrorefiection. At a given energy E' > 0, since 

Ue = {jlTI + E)/{jlTI -E)>1, 



iyh 



^J-TI - E 
hvi 



\ — v'i, cos^ 9i 



(27) 



becomes purely imaginary with Imikyh) < for a certain 
range of 9i . This corresponds to an evanescent Andreev 
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mode. Such modes exist for the ranges < 9i < Oie and 
TT — 9iE < 6i < TT, where 9ie = cos~^(l/z/£;). 

On the TI-2 side we may have either a normahy trans- 
mitted electron with momentum {kx,—kz) or an An- 
dreev transmitted hole with momentum {kx,kzh)- The 
longitudinal momenta for the transmitted electron and 
hole are k^ = ^ {uti + E)'^/ {hv2Y ~ s-iid k^^ = 
■\/(^iT7^^'E)27(^t^^)2^^A!|' respectively. [For E — 0, the 
momentum of the hole will be {kx,kz)]. At a given en- 
ergy E, when z/^ > vi/v2, there exist evanescent Andreev 
modes {Im{kzh) < 0) for the following ranges of the inci- 
dent angle di: < 9i < 92e and tt — 92e < 9i < tt, where 
92E = cos^^[vi / {v2iyE)]- With and being the am- 
plitudes for normal and Andreev transmissions, wc can 
write the wave function in TI-2 as 



*2 



V2 



f y i + COS6I2 ^ 

\/l - cos 02 






-ikz z 



1± 

V2 



( 






Vl — cos ( 



'2h 



\ ^/l + co&92h J 



and 



e'^'^- = {k^ + ikz)/yjkl + kl, 
e*«- = {k.,+ikzh)/^kl + k 



2 

zh ' 



(28) 



function on SC looks like 
1 

*3(2/) = ^ 















< -V" ) 







^{-ikR-ki)y 



where e^^ 



^ {E-^^Al-E^). (30) 



C. Wave EYinction for p^-wave SC and d = z 

For a p-wave SC with f{k) = ky/kF (where fikp = 



\/2rniIsc) and d = z, the dispersion is anisotropic. For a 
given E and 9i on TI-1, the longitudinal momentum on 
the SC side (ifc^ -I- ikj) is given by 



= ±kf 



\ 



where Fi = 



l-T^±i^ 



^0 



^^'sc 



-r2, 



12 — 



^0 



(31) 



and we choose the ± signs so that Im{kysc) > 0. The 
wave function is given by 



B. Wave Function for s-wave SC 



For the case of an s-wave SC, as mentioned earlier the 
pair potential is isotropic and k independent as shown in 
Eq. (5). The momentum along x will be equal to k^ (the 
same as in the TIs) due to the translational invariance 
of the system in the x direction, while the longitudinal 
momentum kysc will be given by 














w 



where w 



U_ 
V2 



^{ikR-ki)y 



I 
1 

1/7* 

V 

/^nkysc/ki. 



E + h^k^/{2m)- fisc 



^{-ikii-ki)y 



(32) 



'iySC 



±kf 



\ 




(29) 



where we choose the ± signs in such a way that 
Im{kysc) > 0. Of the four possible solutions for kysc 
in the above equation, we must choose the two for which 
the wave functions decay as ?/ — >■ 00; they can be writ- 
ten in the form kysc = i^fl + ikj where kji and kj are 
positive. At energies in the gap {\E\ < Aq), the wave 



D. Conductance 

From the different reflection and transmission ampli- 
tudes calculated we can compute the four probabilities 
Rn - IrjvP, Ra = \rA?, Tn = |tw|', and Ta = \tA? as 
fimctions of 9i and E. The boundary condition Eqs. (24) 
and Eqs. (25) imposed conserves both probability and 
charge currents at the junction. The conservation of the 
probability current at the junction implies that 

i;isin6'i = fi(sin6'ii?7v + sin6'i;ji?yi) 

+V2{siTi92TN + &in92hTA) (33) 
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for each value of 9i. From charge current conserva- 
tion we can write down the charge current on SC side 
in terms of the charge currents on TI-1 and TI-2 as 
J3,tot = Ji,m - J2,out- The incoming charge current 
along the y direction on the TI-1 side is equal to Jijn = 
ewi[sin0i(l — Rn) + sin 6*1/1 while the total outgoing 
charge current along the — z direction on the TI-2 side 
is equal to J2.out = ev2[sin6'2r/v — aii'i02hTA]- Using the 
probability conservation Eq. (33), this can be written as- 



(34) 



Essentially, we have written the charge currents on all the 
three sides in terms of only the scattering probabilities 
in TTl and TT2. 

The above discussion needs to be modified if there is 
an evanescent Andreev mode on TI-1. When kyh be- 
comes purely imaginary, the term proportional to va in 
Eq. (26) becomes a decaying wave (rather than a plane 
wave) and therefore does not contribute to the proba- 
bility and charge currents along the y direction. Hence 
Eqs. (33) and (34) will not contain the term proportional 
to Ra- The same argument can be repeated if there is 
an evanescent Andreev mode on TI-2; then Eqs. (33) and 
(34) will not contain the term proportional to Ta- 

To obtain expressions for the various conductances 
from the above currents, we assume that a voltage bias 
V is applied on the TI-1 side maintaining the TI-2 and 
SC at the same potential. Namely, we choose the mid- 
gap energy on the SC side as E = and maintain the 
Fermi energy on TI-2 at zero and the Fermi energy on 
TI-1 at eV. The differential conductance Gi = dli/dV 
is then the derivative of the current measured on side i 
{i — 1,2,3) with respect to V. Integrating the various 
currents over the angle of incidence 9i leads to the follow- 
ing expressions for the differential conductances at zero 
temperature. 



V. CONDUCTANCE RESULTS 

In this section, we compute the sub-gap charge con- 
ductance numerically for a chosen set of parameters for s- 
and p-wave superconductors as worked out in Sec. IV B 
and IV C. The plots of the sub-gap conductance G3 as 
a function of the bias voltage E are shown in Fig. 2 and 
Fig. 3 for the choice xi = X2 = 0. We find that in ac- 
cordance with conventional NM-SC junctions, G3 shows 
peaks at the gap-edge for the s-wave SC and at mid- 
gap for the p-wave SC. These peaks get sharper with 
increasing barrier strengths and for large X3i G3 reflects 
the density of states of the superconductor. Note that 
the sharp mid-gap peak for the p-wave case indicates the 
localized mid-gap edge states. However, in contrast to 
conventional NM-SC junctions, the peak height does not 
reach a value of 2Go for large X3 if Xi = X2 = 0. This be- 
havior is reminiscent of the graphcne-SC junctions and is 
a reflection of the Dirac-like properties of the electrons on 
the TI side^^. To elucidate this property further, we are 
going to turn on xi and X2 and study the behavior of the 
sub-gap conductance as their function. Before resorting 
to such a study, we plot Gi.2,3 for sub-gap voltages and 
for xi = X2 = in Figs. 4 (a-b). As discussed earlier, 
these satisfy the relation Gi = G2 -|- G3. 



0.001 
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-0.5 
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Gi{E) = dei[sm9i{l^RN)+smeihRA], 

Gn r V9 

G2{E) = dOi —[suie2TN-sin92hTA], 

2 Jo VI 

GsiE) ^ GolE / de,[sm0ihRA + -sm92hTA] , 
Jo VI 

, „ 2e2VP>T/ , , E 

where Go = — ; and 7b = 1 H . (35) 

h hvi uti 

Here, we have written Go such that it has the units of 
conductance and we have chosen /iy/ > Ao. The factor 
of ^Ti + E m the product Gq^e comes from the linear 
density of states of the incoming electrons on the TI-1 
side. Finally, for the case V2 = vi and ^ti ^ E (when 
we find that 02, 0ih and 02h are all equal to 6*1), we find 
that the conductances must lie within certain bounds: 
< Gi < 2Go, -Go < G2 < Go and < G3 < 2Go. 



FIG. 2: The sub-gap conductance of an s-wave SC in units of 
Go, for xi = X2 = 0. The Gs-axis is shown on a logarithmic 
scale to make the distinction between different lines clearer. 



Edge States 

To understand the variation of the sub-gap tunneling 
conductance G3 with the barrier potentials xi and X2, 
we first discuss the localized edge states of the super- 
conductor near the line junction. We note that a SC 
with a boundary along y = 0, depending on its pairing 
symmetry, may exhibit sub-gap localized edge states at 
the boundary. These states appear at the edges (mid- 
dle) of the superconducting gap for s- (p-) wave pairing 
symmetry^°. To understand this, let us first consider a 
system consisting of only a SC with a hard wall at y = 0, 
and study the bound states which can occur in that case; 
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FIG. 3: The sub-gap conductance of a p^-wave SC in units of 
Go, for d = £ and Xi = X2 ~ 0. 
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(a) s-wave SC 
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(b) py-wave SC with d = z 

FIG. 4: The differential conductances Gs into the SC, G2 into 
the TI-2 and Gi back to the TI-1 in units of Go, for (a) an 
s-wave SC, and (b) a py-wave SC with d = z. 



later, we shall discuss the effect that such bound states 
have on the conductance when the SC is coupled to the 
two TIs. 

The wave function of a bound state will have a factor 
^i{k^x+kyy~Et/h) ^ where and E are real and Im{ky) > 
0. Given some real values of E and k^, dispersion rela- 
tion will generally have four solutions for ky of the form 
k, k*, —k and —k*. Of these we choose the two solu- 
tions whose imaginary part is positive. Let us call those 



solutions k± = ±kf! + ikj^ where k^ and fc/ are posi- 
tive and denote the corresponding four-component eigen- 
states of the Hamiltonian by and "I*- for a particular 
component of the spin. If these eigenstates are identi- 
cal, then we can take their superposition with opposite 
signs to obtain a wave function which is proportional 
to ^+sm{kjiy)e^^''^^~^*/^^~'''y; this vanishes at ?/ = 
thereby satisfying the hard wall boundary condition and 
giving us a bound state. In general, for a given value of 
E, there will be a particular value of kx for which there 
are two eigenstates which are identical and can therefore 
be superposed to give a bound state; hence the bound 
states will have a dispersion in which E' is a function of 
kx- Further, there will be two such states for a given E 
and kx corresponding to the spin degree of freedom. 

Now we discuss the effect that the existence of bound 
states on the SC-side may have an effect on the conduc- 
tance when the SC is coupled to two TIs. In general, an 
electron incident on TI-1 penetrates into the SC up to 
some distance (cx l/ku) and the amplitude (hence the 
transmission probability) for such a penetration can be 
suppressed by increasing the barrier strength X3. How- 
ever, if an electron is incident with an energy which is 
exactly equal to the bound state energy on the SC side, 
there will be a resonance in the transmission to the SC 
side and the barrier xs does not affect the transmission. 
A flat dispersion for the bound state spectrum therefore 
means that the conductance should be enhanced when 
the bias matches exactly the bound state energy on SC 
side. We shall soon see that the bound state spectrum is 
flat for certain cases. For larger barrier strength X3, the 
peak is more visible since a larger barrier will suppress 
transmission at energies other than the bound state en- 
ergy. We find such peaks in our numerical calculations 
for (i) an s-wave SC at E = ±A (see Fig. 2) and (ii) a 
p-wave SC at = (see Fig. 3) as expected^", 
s-wave SC: The wave function ^3 for an s-wave SC is 
given by Eq. (30). We find that the wave function ^3 
consistent with the hard wall boundary exists only when 
E = ±A. The dispersion for the bound states is flat 
E — ±A for any kx- These bound states are responsible 
for the sub-gap conductance peaks at i? = ±A in Fig. 2. 
At i5 = A, 77 = and hence ^3 and can be written 

as a linear combination of the spinors [10 1]'^ and 
[0 1 - 1 0]^. So if we multiply Eqs. (24) from the left 
by the orthogonal row vectors [1 0—1] and [0 1 10], 
the boundary condition reduces to four equations with 
four unknowns tjy, tA, rrf and which can easily be 
solved. (The details of a similar calculation for a p-wave 
SC are presented in the Appendix). The important point 
to note here is that the barrier strength xa drops out of 
the problem. Let us choose Xi = X2 = for simplicity. 
We find that iyv = = and 



\rA\' = 



1 — cos ( 



1 



l/l) 



(36) 



A similar calculation yields the same expression for the 
scattering probabilities for the case E = —A. However 
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the difference to note is that 9ih is different at energies 
E = ±A since cos^i^^ = cos 6*1 {fixi + E)/{^ti — E). 
Numerically, we find that the conductance G3 approxi- 
mately obeys G3(A) = G3(— A) in the limit use ^ ^■ 
p-wave SC: Let us study a p-wave SC with the choice 
f{k) — ky/kp and d ^ z. For this case, the wave function 
^3 given by Eq. (32) is consistent with the hard wall 
boundary only when w = —w* which simplifies to 



E 



h\kl 



2m 







(37) 



We find that the bound state dispersion is flat: E ^ 
for all kx- [However when /(fc) oc ky + iak^, with a 
non-zero a, the bound states have a dispersion owing 
to TRS breaking]. From the calculation detailed in the 
Appendix, we see that for E~0,rA=tN^O, 



1- COS 26*1 

1 + sin^ 01 + cos2 6*1 sin [2(xi - X2)] ' 

and |7'Arp = 1 — Remarkably, ji^ip and hence the 

mid-gap conductance G3(0) depend only on the differ- 
ence Xi ~ X2- [This is in contrast to the conductance of 
a junction of just two topological insulators (i.e., our set- 
up without the SC side) where the total barrier strength 
seen by the electrons is X1+X2; then the conductance will 
only depend on xi +X2]- The important result that 
depends only on the combination xi ~ X2 (and hence G3 
also, since ta = 0) at zero bias can be understood quali- 
tatively as follows. An electron incident on the junction 
from the TTl side sees a barrier strength of xi on that 
side. It then transmits to the TT2 side with an ampli- 
tude tA as a hole; such a hole sees a barrier strength of 
— X2 on that side, since the potentials seen by electrons 
and holes have opposite signs. Hence the total barrier 
strength seen is given by xi ~ X2\ hence this is the pa- 
rameter which appears in Eq. (38). 

The numerical results shown in Fig. 5 highlight these 
features. To illustrate the behavior of conductance as a 
function of different barrier strengths, we choose xi = 
X2 = X3 = X- At non-zero bias, the conductance has 
periodic oscillations with a decaying envelope as shown 
in Fig. 5 (a). The periodic oscillations are due to the 
barriers xi and X2 on the TI surfaces while the decaying 
envelope is due to the barrier xs on the SC side. At 
zero bias (i.e., = 0), the barrier on the SC side no 
more plays a role in determining the conductance as we 
saw analytically earlier. However, the barriers xi and X2 
on the TI side together determine the value of the ZBP. 
The numerical results in Fig. 5 (a) for X2 ~ Xi and (b) 
for X2 = — Xi show that the value of G3 at E' = is a 
function of Xi ^ X2 only. 

For a p-wavc SC, when d = x ot d = y, we find that 
the dependence of the conductance on the bias is different 
from that for d = z shown in Fig. 3. For the case d ~ y, 
we again find that r^i = t^r = while 



\tA\ 



1 - COS 26*1 



1 + sin^ 9i - cos2 01 sin [2(xi - X2)] 



(39) 




(a) Xi 



X2, 




(b) XI 



X X-i /Co Xq 



-X2, E = 



FIG. 5: (a) Conductance of a py-w&ve SC with d = z at 
different energies in the SC gap. (b) Conductance of a py- 
wave SC at is = for different spin pairings: d — z (blue 
dashed hne), d = x (green solid line) and d = y (red dot- 
dashed line). Note that G3/G0 reaches 2 a.t x ~ mr /2 + 3tt /8 
for d — z and at x = ?^7^/2 + 7r/8 for d = y, as can be shown 
from Eqs. (38) and (39). 



(The qualitative reason for the dependence on xi ~ X2 is 
similar to the one given above for d ^ z). For the case 
d ~ X, wc find that ^ tpf = 0, while 



\rA\^ = sin^ 01 



(40) 



and {tnI'^ = cos^ 0i are independent of xi and X2- This 
independence stems from following reason: an electron 
incident on the junction from the TI-1 side sees a barrier 
strength of xi- When it is Andreev reflected as a hole 
back to that side, it sees a barrier strength of — Xi- Hence 
the total barrier strength that is seen is zero, regardless 
of the value of xi- Further, there is no dependence on X2 
since the electrons does not transmit to the TI-2 side at 
all; both Ia and In are zero. 

These analytical results are confirmed by our numeri- 
cal calculations. This is illustrated in Fig. 5 (b) where G3 
is plotted as a function of xi — X2 at i? = for the three 
different forms of d. Note that for d~y and z, G3 varies 
with a period of 7r/2 in the variable x — (xi ~ X2)/2 in 
Fig. 5 (b). This is in agreement with Eqs. (38) and (39) 
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which are periodic in Xi ~ X2 = i"- We note that our re- 
sults demonstrate that the barrier dependence of the sub- 
gap tunnehng conductance of junctions of TI and triplet 
superconductors can be used to determine the direction 
of the direction of d, or equivalently, the direction of the 
spin of the Cooper pair. This property is in complete 
contrast to the conventional NM-SC junctions and stems 
from the spin-momentum locking of the Dirac electrons 
which breaks the rotational symmetry in the spin space 
necessary to distinguish between the different orientation 
ofd. 

For a p-wave SC, we can in general take f{k) oc ky + 
iakr^. The choice a = corresponds to Py-SC and the 
bound state dispersion is flat as discussed earlier. When 
a is changed slightly from zero, we find that the bound 
state dispersion develops a small non-zero slope. Now, 
the ZBP which existed for a = broadens and the peak 
value also decreases; the peak eventually disappears as a 
is increased further. 



VI. SPIN CURRENTS 

In this section we will study the spin currents on the 
SC side for s- and p-wave cases. For a = x,y, z, the a 
component of the spin density is given by 



Pa 



2 3 



(41) 



in the four-component language. The operator ap- 
pears in Eq. (41) because a hole with spin-up (down) 
corresponds to a missing electron with spin-down (up). 
The spin current Ja corresponding to the spin density 
Pa can be calculated on the SC side using the equations 
of motion Eq. (11) and (12) and the continuity equation 
dtPa+^ ■ Ja = 0. The equation of continuity implies that 
in a steady state, the component y ■ Ja is independent of 
the y coordinate and can therefore be evaluated at any 
convenient point like y = or 00. We are interested in 
only the y component of the a-spin currents since only 
this component contributes to the spin density that is 
injected into the SC. The differential spin conductance 
Ga corresponding to the spin density pa is obtained by 
integrating over the spin current contributions y ■ Ja from 
all the angles of incidence < 61 < tt : 



where Gn 



hv'j 



PSC Jo 



dOi y ■ Ja 



(42) 



Here Gn has been defined in 



such a way that the term in the brackets in the above 
equation becomes dimensionless. 

In a normal metal, i.e., setting Aq = in Eq. (11), we 
find that 



Ja 



2m 



(43) 



In a SC, the spin current Ja is a sum of two parts: the 
part Ja,N which is independent of Aq given by Eq. (43), 
and the part Ja,A proportional to Aq. Ja,N is the ex- 
pression for the a-spin current in a normal metal carried 
by either electrons or holes; in SC it gets contribution 
from the BdG quasiparticlcs. For an s-wave SC, we find 
that the x- and z- spin currents (i.e., Ja for a — x,z) 
are entirely given by Eq. (43) and do not contain the SC 
pair potential Aq {Jx/z,a — 0). Since the wave function 
^'3 decays exponentially as y — >■ 00, the spin current go- 
ing into the SC is zero at y ^ 00 (and therefore at any 
value of ?/) if a = x or z. But the y-spin current in ad- 
dition to the expression in Eq. (43), also contains a part 
proportional to Aq given by- 



y ■ Jy-A 



dy' Re[^l{x,y')T-^3{x,y')]. (44) 



The derivation of this is similar to that of Eq. (18). Hence 
the total spin current is given by the sum of Eqs. (43) and 
(44), and this need not vanish. It is easiest to evaluate 
this at y = where Eq. (44) vanishes. For a = y, after 
tedious calculation, wc find that Eq. (43) simplifies to 



y ■ Jy 



h^kE 
2m 



{1 + \w\^) Im[tlt2 - t;t4]. (45) 



The above quantity is non-zero for the scattering of an 
electron incident at energy E, at an angle 61 on TI-1. 
But the y ■ Jy currents at incident angles 9i and -k — 9i 
add up to zero, thus contributing nothing to the differ- 
ential spin conductance Gy. To summarize, all the spin 
conductances G^ are zero for the s-wave case. 
For a p-wavc SC, the a-spin current is given by 



V • J„ A = Aq Re 



a - [V'V''(cr-CT)icr^V0 
-Vi:\d-a)iaya''(j)]\. (46) 



We will now consider all the nine different possibilities 
corresponding to d = x,y,z and a = x,y, z. Table I sum- 
marizes the results for different spin currents, while Ta- 
ble II summarizes the spin conductance results (obtained 
by integrating out spin currents over all incident angles 
Oi) for different cases. For the four cases corresponding 
to (c?, a) = {x,x), (x,y), {z,y) and [z,z), we find that 
V • Ja,A is a total derivative, which implies that Ja.A has 
a local expression in terms of \l/3 (analogous to the last 
terms in Eq. (19)). Then one can evaluate both Ja.A and 
Eq. (43) at y = 00; they vanish there because goes to 
zero exponentially. Hence no spin current enters the SC 
in these four cases. In the remaining five cases, we find 
that V • Ja.A is not a total derivative; hence J a, a has an 
integral expression (analogous to Eq. (44)). Since this 
vanishes at y = 0, one can find the total spin current by 
evaluating only Eq. (43) at y = 0. 

d ^ X : As argued earlier, no a-spin current correspond- 
ing to a = a;, y enters the SC. Jz on the other hand has 
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p-wave 


p-wave 
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Eq. (48) 


Eq. (48) 


Jy 


Eq. (45) 





Eq. (49) 
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Eq. (47) 


Eq. (47) 






TABLE I: Expressions for spin currents for s-wave and p-wave 
pairing for a particular angle of incidence. 
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p-wave 


p-wave 


p-wave 
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non-zero 


non-zero 
















Gt 















TABLE II: Spin conductances Ga for s-wave and py-wave SC. 



an expression given by 



y-Jz 



Am 



■(l + IH^) [\tl\'~\t2\'-\t3\' + \h\']. (47) 



This is the expression for z-spin current for electron in- 
cident at a given energy E and angle 9i on TI-1. The 
z-spin conductance on SC is obtained by summing up the 
currents at all possible incident angles < di < tt. The 
z-spin conductance Gl turns out to be zero due to the 
exact cancellation of the y ■ Jz at angles Oi and tt — . 
d = y : All the three spin currents could have a non-zero 
value in this case (see Table I). The expression for z-spin 
current in this case is given by Eq. (47) and the total z- 
spin current at any given energy E in the gap goes to zero 
by the same logic that follows Eq. (47). The expressions 
for the x/y-spin currents are 



y ■ Jx 

y- Jy 



2m 



2m 



{1 + \w\'') Re[tlt2 - t;t4], (48) 
■{l + \w\^) Im[tlt2-t;u]. (49) 



Integrating the spin current over 0i, we find that the y- 
spin conductance is zero while the x-spin conductance is 
non-zero and has features similar to the d — z case as 
discussed below. 

d = f : In this case as argued earlier, y/z-spin currents 
are zero. The expression for x-spin current is given by 
Eq. (48). The spin conductance which is integrated over 
a;-spin current is non-zero and has some feature as shown 
in Fig. 6 for a typical case. 

If we compare Tables I and II, the non-zero ?/(z)-spin 
current for an electron incident at angle 6i in Table I 
adds up with the y(2)-spin current for an electron in- 
cident at angle i: — 9i to give zero in the total y{z)- 
spin conductance. This can be seen from the symme- 
try of the equations of motion and the boundary con- 
ditions under the transformation x — —x. For both 



TI-1 and TI-2, the equations of motion (Eq. (10)) re- 
main invariant under the transformations: — > cr^^'i 
and 'f i T^a^'i'i accompanied by a; ~x. Both the 
boundary conditions Eqs. (24) and Eqs. (25) are invari- 
ant under these two transformations. However, for s- 
wave SC and the p-wave SC with d = x the equations of 
motion (Eq. (11) and Eq. (12)) are invariant under the 
transformation 'i's — > r^cr^^'s. And for p-wave SC with 
d ~ y, z the equations of motion (Eq. (12)) are invari- 
ant under the transformation ^'3 — > cr^^'a. Now, we can 
easily see from Eq. (43) that under both the transforma- 
tions {'i'i — > (T^'i'i and '^i r^cr^^'i), y ■ is invariant 
while y ■ Jy and y ■ Jz change sign. Since x — >■ —x is same 
as 01 — TT — 01 , this explains the exact cancellation of y- 
and 2-spin currents for di and tt — 0i . 

Wc can relate the spin currents which are non-zero in 
Table II to the physical quantities on the other side of 
the junction, i.e., on the TI-1 and TI-2 sides. Using the 
boundary conditions in Eqs. (25), we can rewrite the x- 
spin current at the junction as 



y 



J 



'SC 



-[uiflr'^i 



(50) 



for a pj,-wave SC with d — y or z. We thus see that the 
x-spin current on the SC side of the junction is linearly 
related to the steady state charge densities on the TI-1 
and TI-2 sides of the junction evaluated at the junction. 



Spin Conductance 

The x-spin conductance shows an unusual satel- 
lite peak (SP) in addition to the ZBP for the cases d = y 
and d = z when f{k) = ky/kp (see Fig. 6). This SP 
merges with the ZBP for X3 }^ 20. The SP is observ- 
able only if there is an appreciable penetration of a plane 
wave state with non-zero energy from the TI-1 into the 
SC, and this can occur only if xs is not too large. The 
location and height of the SP change with xi and xi in 
a periodic way. To highlight these features, we choose 
X\—X2 — Xz = X S'lid show the spin conductance G% as 
a function of the bias and the barrier strength x in Fig. 7 
as a contour plot. 

Wc find, using Eq. (50), that the x-spin current in the 
superconductor can be written in terms of the charge 
densities in TI-1 and TI-2. In the case V\ ~ V2, we find 
that the spin current is proportional to \I'|r^\E'i — 
which is given by 



= 1-1- RjM + RN,Ph — JiA^lA — Tn 



TaT2aA^^) 



where 

RN,Ph 



cos^ 9i — ve cos 9i \lv% cos^ 6*1 — 1, 
\i'ECOsei\, 



(52) 
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in the range of {E,9i) for which there exist evanescent 
Andreev modes on TI-1 and TI-2, i.e., when t'£;cos0i > 
1. When the Andreev states are plane wave states on TI- 
1 and TI-2, Via and T2A both are equal to 1 in Eq. (51). 
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FIG. 6: Spin conductance for a Py-wave SC with d = z and 
Xi ~ ~ 0- There is a satellite peak at a positive value of 
the bias which merges with the ZBP in the limit of large X3- 



X 




FIG. 7: Spin conductance for a py-wave SC with d — z and 
Xi = X2 = X3 = X- Both the location and the height of the 
SP oscillate and decay with increasing x- 

Using the identity in Eq. (51), we find the contributions 
to the spin conductance from the different scattering 
amplitudes on the TI-1 and TI-2 sides. It is interesting to 
note that in contrast to the expressions in Eq. (35) for the 
charge conductances where only the scattering probabil- 
ities in the TI-1 and TI-2 appear, the spin conductance 
on SC side depends both on the phase and the ampli- 
tude of the reflection amplitude rjv- We show in Fig. 8 
the contributions to the spin conductance from the dif- 
ferent terms in Eq. (51). From the figure, it is evident 
that the Andreev scattering term {Ta^2A — Ra^ia) and 
the term Rn.pii (which is proportional to rjv) contribute 
the most to the SP. Further, to understand the origin of 



the SP, we look at the spin currents due to incident elec- 
trons with different energies and angles of incidence 0i. 
In Fig. 9, we show spin conductances at different E and 9i 
as a contour plot. This contour plot shows that the ZBP 
and the SP in the spin conductance get contributions 
from different angles 61. In contrast to the ZBP in the 
spin conductance which gets contributions mainly around 
normal incidence (^i ^ 7r/2), the SP gets major contri- 
butions from electrons incident at glancing angles {9i ^ 
and 6*1 ~ tt). This is highlighted in the contour plot by 
ellipses. The green (dot-dash) lines, 61 = cos~^(±l/i^£;), 
separate the evanescent Andreev modes on TI-1 and TI-2 
from the Andreev plane wave states. The largest contri- 
bution to the SP comes from the range of 9i for which 
there are evanescent Andreev modes; this is shown by 
the pink (dark) region in the contour plot. In the same 
region, we find that the phase of r^y changes rapidly from 
— TT to TT (going through 0) as E is varied at a fixed value 
of 9i. Such a feature in the phase of shows up as a 
peak in the term RN,Ph in Eq. (51) when E is varied at 
a fixed 9i. Further, the SP in spin conductance merges 
with the ZBP in the limit /isc ^ 00. 




FIG. 8: Spin conductance and the contributions from different 
terms in Eq. (51) for a p^-wave SC with d = z and xi = 
X2 = X3 = 0. The phase term Rn^ph and the Andreev terms 
TAr2A - RaTia contribute to the SP. 



There is a striking asymmetry about E ^ in the 
spin conductance in Fig. 6; this mainly arises because 
the evanescent Andreev modes can exist only for i? > 
and are absent for E < 0. There is also an asymmetry 
due to a finite difference in the density of states of the 
incident electrons at energies E and —E (see Eq. (35)) 
and the lack of invariance of Eq. (32) under E —E. A 
small asymmetry is also visible in the charge conductance 
shown in Fig. 3. 
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FIG. 9: Contour plot of spin current y ■ (in units of Jq — 
h'^ kp /{2m)) as a function of the energy and incident angle of 
the electron on the TI-1 side, for a p^-wave SC with d — z 
and Xi = X2 = X3 = 0. The dominant contribution to the 
SP comes from the pink (dark) regions which are bounded by 
the green dot-dash lines. See text for details. 



VII. SUMMARY AND DISCUSSION 

In this paper, we have considered a junction between 
two TIs and a SC which can be either s-wave or p-wave; a 
p-wave SC may have different spin pairings and we have 
studied each case. We have formulated the most general 
time-reversal invariant boundary condition at the junc- 
tion; this consists of three barrier parameters. We have 
studied the dependence of the differential charge and spin 
conductances on the barrier parameters, the ratio of the 
TI and SC chemical potentials to the superconducting 
pair potential, and the voltage bias. 

Our main results can be summarized as follows. 

(i) In accordance with conventional NM-SC junctions, 
we find that the charge conductance G3 shows peaks at 
the edges of the SC gap for s-wave SC; for a p-wave SC, 
there is a peak at the zero bias. However, in contrast 
to conventional junctions, the heights of these peaks do 
not reach the conventional value of 2Go, if Xi = X2 = 
0; this feature is a signature of the Dirac nature of the 
TI quasiparticles. Further, the height of these peaks is 
independent of X3- The conductance peaks arise due to 
the presence of bound states on the SC side; as a result, 
an electron coming in from the TI side with the same 
energy as the bound state can resonantly enter the SC, 
and then convert into Cooper pairs deep inside the SC. 
(The phenomenon of a conductance peak resembles the 
transmission resonance which occurs in other systems; 
see Ref. 41 for a few examples). 

(ii) For a p-wave SC, we studied the conductance as a 
function of barrier strengths. At non-zero energies the 
conductance both oscillates and decreases with increasing 
barrier strength. This demonstrates both the Dirac and 
the Schrodinger nature of the electrons in different parts 
of our set-up in the following sense. The transmission of 
Dirac electrons oscillates with the barrier strengths in the 



TI-1 and TI-2 (as can be seen by the presence of cosxi 
and sinxi terms in the boundary conditions (24) and 
(25), for i = 1,2), while the transmission of Schrodinger 
electrons through the barrier in the SC decays exponen- 
tially with the barrier strength xs- This indicates that 
for Xi ^ X2 = the charge conductance in these junc- 
tions docs not reach the maximum value of 2Go which is 
found in conventional junctions. 

(iii) For a p-wave SC, the dependence of the conductance 
on the barrier parameters (specifically on the difference 
of the parameters in the two TIs. xi ~ X2j Sut zero bias) 
varies with the direction of d, i.e., the spin pairing of the 
Cooper pairs. This is due to the spin- momentum locking 
of the electrons in the TIs and the fact that the geomet- 
rical arrangement of the two TIs completely breaks the 
rotational symmetry. This is quite different from what 
happens in a junction of normal metals with a p-wave 
SC; since a metal is isotropic in spin space, the conduc- 
tance into a SC would not depend on the direction of d. 
Thus these junctions are expected to provide a test bed 
for determining the orientation of d and hence the Cop- 
per pair spin direction of triplet SC without application 
of external magnetic fields''^ . 

(iv) The spin conductance into the SC for different pair 
potentials exhibits several unconventional features. We 
find that there is a non-zero conductance for the x com- 
ponent of the spin only for a p^-wave SC in which d = y or 
z. (The fact that all the other spin conductances are zero 
can be shown using a parity symmetry which is present 
in our system). The non-zero spin conductance arises 
because the TIs break the spin rotation symmetry, and 
it docs not occur for a NM-SC junction. Most remark- 
ably, the spin conductance shows a satellite peak away 
from zero energy. We have provided an explanation of 
this satellite peak by relating the spin conductance on 
the SC side to the charge densities on the TI sides. The 
charge density turns out to depend on the phase of the 
normal reflection amplitude which shows a peak at non- 
zero bias. The satellite peak in the spin conductance is 
a distinctive feature of our system and it would be inter- 
esting to look for this experimentally. 

The simplest experimental verification of our theory 
would be to measure the maxima of the charge conduc- 
tance peaks for both s- and p-wave SCs. A prediction 
of our theory is that the values of these peaks reach 2Go 
only for some non-zero barrier strengths. It will also be 
interesting to measure the barrier potential dependence 
of G3 . The barrier potentials can be tuned by applying 
a gate voltage along the junction in our system. It may 
be difficult to tune xi; X2 and X3 separately. Looking 
at Fig. 1, however, it is apparent that if a gate voltage 
is applied from above the junction, it will have a larger 
effect on xi and X3 than on X2 (since it is further away 
from the TI-2); hence changing such a gate voltage will 
change the value of xi ~X2- Hence our important predic- 
tion that the zero bias peak depends on xi — X2 for d ~ y 
and z but does not depend at all on any of the Xi for 
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d = X is something which can be tested by changing the 
gate voltage. A consequence of this would be that mea- 
surements of G3 can be used to determine the direction 
of d for a given p-wave SC by attaching it to a junction 
of two TIs. 

Recently, there has been much excitement about 
a zero-bias peak (ZBP) observed in a number of 
experiments in semiconducting/superconducting 
nanowires''^"^^ . This peak is believed to be due to 
a Majorana fermion mode, and it has a close connection 
to the mid-gap (zero energy) states which are present 
at the boundary of a p-wave SC^"*""^^. The last few 
years have witnessed intense theoretical activity in the 
area of Majorana modes at the ends of one-dimensional 
systems^^"^^. In our system, we have studied the edge 
states which occur at £^ = in a p-wave SC with a 
hard wall. However, although these states give rise to a 
zero bias peak in the charge conductance, they cannot 
be called Majorana modes for the following reason. 
A Majorana mode must have a real wave function so 
that it remains invariant under complex conjugation. 
However, our zero energy states carry finite momentum 
kx and do not have real wave functions. Under complex 
conjugation, a state with momentum kx changes to 
a state with momentum —kx and therefore does not 
remain invariant. 

Finally, we point out that in this work we have studied 
the charge and spin conductances of a junction between 
two TI surfaces and one SC. It is natural to ask what 
would happen if there was a junction of a single TI sur- 
face and an SC. A peculiarity of this problem would be 
that a single TI has only one spin degree of freedom for a 
given value of the energy and momentum, while a SC (or 
metal) has two spin degrees of freedom. Hence boundary 
conditions such as the ones given in Eqs. (24-25) would 
generally be inconsistent since they provide eight equa- 
tions for six amplitudes (normal and Andreev reflection 
in the TI and four amplitudes on the SC side). For our 
system with two TIs and a SC, there is no such mis- 
match. This enables us to find all the amplitudes for 
all values of the system parameters, and the amplitudes 
always satisfy the conservation of probability. 
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Appendix 

Here, we shall present a calculation to show how dif- 
ferent scattering amplitudes become independent of the 
barrier X3 on the SC side when the SC side hosts bound 
states. In particular, we choose to discuss the case of 
Py-SC with d = z, and vi = V2 in the TIs. 



From the Hamiltonian in Eq. (4), it is easy to see that 
for a Py-SC with d = z, w defined in Eq. (32) satisfies 
the 2x2 matrix equation 

LV 2m / kp i \w J \ ^ / 

(53) 

At i? = we can see that the two rows of this equa- 
tion imply that w'^ = —1. Further, ky ~ kysc given by 
Eq. (31) can be written as ky ~ ka + iki with kj > 0. 
With this, it is easy to see that w = —i. Then, Eq. (32) 
shows that the wave function on the SC side ^3(2/) and 
its derivative dy'^sd/) both reduce to a linear combina- 
tion of only the two spinors vj/^^ = [1 — i]'^ and 
^'pj, = [0 1 — i 0]"^. Also, the third term on the left hand 
side (LHS) of the second boundary condition in Eqs. (25) 
will be a linear combination of the spinors "^Ip-^ and '^p^ 
for a Py-SC with d = z. Hence the LHS of both the 
equations in the boundary condition in Eqs. (25) arc lin- 
ear combinations of the spinors ^p-j- and '^pi- Since both 
^'Jl and vl'jj^ are orthogonal to ^p-|- and 'i'pi, multiplying 
by "iip^ and from the left in Eqs. (25) makes the LHS 
equal to zero; hence the equations become independent of 
X3. We then get the following four equations containing 
four unknowns r^r, rA, and Ia, 

^l^[Mixi)^i+M{-X2)^2] = 0, 

*JjM(xi)*i+M(-X2)*2] = 0, 

*J^a-r^[M(Xi)*i-M(-X2)*2] = 0, 

*^^a^T^[Af(xi)*i-A/(-X2)*2] = 0, (54) 

where ^1 and 5*2 are evaluated at y = 0. When 'i'l and 
\I'2 from Eqs. (26) and (28) arc substituted in the above 
equations, we get four equations in four unknowns r^, 
rA, tN and Ia which can be cast in the matrix form: 
M-X = A where 

9(Xi) -w^eif'(Xi) HX2) wh*(x2) 

jyj ^ -7eiS*{-Xi) wg{-xi) + X2) -™/i(f - X2) 

-7eiS*(-Xi) -wg{-Xi) j/i(f+X2) -jw/i(f-X2) 

9(Xi) ^'YnQ'iXi) -HX2) wh*{x2) 

X = [tm rA tN Ia]^ , 

A = [-g*(xi) -ia*{^ + xi) -i9*{^ + xi) -g'ixi)] , 
and g{x) = cos x + 17^1 sin x, je, ^ ~ie''^\ 
and h{x) ~ ^1 + cos 61 cos x + « yl— cos^ sin x. (55) 

This set of linear equations yields rA = = 0, 



1 + sm t'l+cos^^i sm[2(xi-X2)J 

and I^Arp = 1 — \tA\'^- The cases d = x and y follow a 
similar calculation, and the results are given in Eqs. (39) 
and (40). 

For a s-wave SC, a similar calculation can be done at 
E = ±Ao; once again, the existence of a bound state 
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on the SC side with a hard wall boundary makes two eigenspinors and one can use them to reduce the problem 
eigenspinors of the SC Hamiltonian identical. Therefore, to four equations in four unknowns as shown above, 
there arc two other spinors which are orthogonal to these 
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